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We propose a mechanism to evade the Lyth bound in models of inflation. We minimally extend 
the conventional single-field inflation model in general relativity (GR) to a theory with non-vanishing 
graviton mass in the very early universe. The modification primarily affects the tensor perturbation, 
while the scalar and vector perturbations are the same as the ones in GR with a single scalar field 
at least at the level of linear perturbation theory. During the reheating stage, the graviton mass 
oscillates coherently and leads to resonant amplification of the primordial tensor perturbation. After 
reheating the graviton mass vanishes and we recover GR. 


Introduction Inflation [ij is the leading paradigm 
of very early universe cosmology, but its physical origin 
is still mysterious. The generation of primordial gravita¬ 
tional waves is a generic prediction of the inflationary uni¬ 
verse. It leads to B mode polarization in the CMB, and 
provides an important window to the physics of very early 
universe. It was reported that the primordial tensor-to- 
scalar ratio is r < 0.11(95% CL), based on Planck full 
sky survey Q . Several next-generation satellite missions 
(CMBPol, COrE and LiteBlRD) as well as the ground 
based experiments (AdvACT, CLASS, Keck/B1CEP3, 
Simons Array, SPT-3G) and balloons (EBEX, Spider), 
are aimed at measuring primordial gravitational waves 
down to r ^ 10“^. See Ref. 0 for a recent updated 
forecast on these future experiments. 

According to the Lyth bound Q, the tensor-to-scalar 
ratio is propotional to the variation of the inflaton field 
during inflation, i.e. A(j)/Mp ~ f dNy/r/8. The thresh¬ 
old A(f) = Mp then corresponds to r = 2 x 10“^, assumed 
that tensor power spectrum is nearly scale-invariant. The 
sizeable amplitude of the primordial gravitational waves 
requires a super-Planckian excursion of the inflaton, i.e. 
Afj) > Mp. 

In quantum field theory, the naturalness principle tells 
us that the variation of a field cj) over the distance greater 
than the cutoff scale is generally regarded as being out of 
the validity of the theory. In a gravitational system, we 
take the Planck mass as the UV cutoff scale, because 
gravity strongly couples to the matter sector and the 
graviton-graviton scattering violates unitarity above this 
scale. Thus the inflationary prediction may not be reli¬ 
able in the case of a super-Planckian excursion. There¬ 
fore, the detection of the primordial tensor perturba¬ 
tion with its amplitude larger than the threshold value 
r = 2x10“^ has a profound impact on our understanding 
of fundamental physics. It implies that either quantum 
field theory or gravity may be modified in the very early 
universe. 

In this letter, by means of modifying gravity, we pro¬ 
pose a new mechanism to evade the Lyth bound. We 
consider a minimal extension of GR with a non-vanishing 
graviton mass term in the very early universe. Specifi¬ 
cally we propose a model in which the graviton mass is 


proportional to the inflation during reheating. Then the 
coherent oscillation of the inflaton induces that of gravi¬ 
ton mass and gives rise to resonant amplification of the 
primordial tensor perturbation. This is a broad paramet¬ 
ric resonance which includes all long wavelength modes, 
given the graviton mass is much greater than the Hubble 
constant during reheating. After reheating, the graviton 
mass vanishes as the inflaton decays and we recover GR. 

A massive gravity theory The theoretical and 
observational consistency of massive gravity has been 
a longstanding problem, the pioneering attempt could 
be traced back to Fierz and Pauli’s work in 1939 Q. 
However, Fierz-Pauli’s theory and its non-linear comple¬ 
tion, the so called dRGT massive gravity Q, suffer from 
many pathologies Ml- The origin of these pathologies 
is probably the Poincare symmetry of the Stuckelberg 
scalar field configuration. 

Away from the Poincare symmetry, a broad class of 
massive gravity theories have been discussed in the liter¬ 
ature [Imi]. In this letter, we consider a massive gravity 
theory with the internal symmetry [l^ [l^ 

+ ( 1 ) 

where A* is the SO{3) rotational operator, S® are 
three arbitrary functions of their argument, and 
are four Stuckelberg scalars with non-trivial VEVs, 

p° = f{t), p^=x\ z = l,2,3. (2) 

These nontrivial VEVs give a non-vanishing graviton 
mass. At the first derivative level, there are two com¬ 
binations of the Stuckelberg fields that respect this sym¬ 
metry. 

The graviton mass term could be written as a generic 
scalar function of the above two ingredients. 

Due to the internal symmetry p'’ (v?®), there 

are only 3 dynamical degrees of freedom (DOF) in our 
theory, i.e. 2 tensor modes, and 1 scalar mode. In the lan¬ 
guage of ADM formalism or the (3 -|- l)-decomposition of 
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space-time, we find that these two ingredients in Eq. © 
are free from the shift fV® and thus the associated Hamil¬ 
tonian of gravity is linear in fVh This implies that 3 
momentum constraints and the associated secondary con¬ 
straints eliminate 3 DOF in hij, and the number of resid¬ 
ual DOF is thus 3 [l5| . 

Now we apply this massive gravity theory to the early 
universe. To minimize our model, we identify the time¬ 
like Stuckelberg scalar with the inflaton scalar field (p, i.e 
gp = <p. By doing this, we achieve a minimal model of 
massive gravity, in which only the tensor modes receive a 
modification, while the scalar and vector modes remain 
the same as the ones in the single scalar model in GR. 

To be specific, we consider the following action with 
enhanced global symmetry —>• constant • , 


S = 




M2 1 

-^n--g^''d^pd.cp-v{<p) 


-oMlml iP) 




( 4 ) 


where a, P, ip and E are scalar, Si and Fi are vector, and 
'jij is tensor. The vector modes satisfy the transverse 
condition, diS^ = diF^ = 0, and the tensor modes satisfy 
the transverse and traceless condition, y- = = 0. 

Tensor perturbation The action for the tensor per¬ 
turbation reads 



dtdPxa? 





.(9) 


We see that the graviton receives a mass correction. We 
quantize the tensor mode as 

liA^) = X] [afceij(fc,s) 7 fce*'''“' -|- h.c] , (10) 

s=±“' 


where ak is the annihilation operator and eij{k, s) is the 
transverse and traceless polarization tensor which we nor¬ 
malize as 


eij{k,s)F^{k,s') = S^s' ■ ( 11 ) 


where V{tp) is the inflaton potential, the numerical fac¬ 
tor 9/8 is inserted for later convinience, and 5ZS is a 
traceless tensor dehned by 0 


^Zij = 


Ztkzkj 


( 5 ) 


where Z^ is defined by Eq. m with replaced by p, 
Z = Z^^dij, and the summation over repeated indices 
is understood. Noted that the 2nd line of Eq. O is 
the graviton mass term, which does not contribute to 
the background energy momentum tensor. Its non-trivial 
contribution starts from the quadratic action in pertur¬ 
bations. 

As for the mass parameter rri^ (</)), we assume the fol¬ 
lowing scalar field dependence: 


mliP) 


A^2 

1 -I- {p /())*)^ ’ 


( 6 ) 


where </* is the inflaton field value at the end of inflation. 
Without loss of generality, we assume p = 0 is the mini¬ 
mum of the potential at which the inflaton settles down 
after reheating. 

As usual, we consider a flat FLRW background, 


= —dt^ -I- aPdx^. (7) 


Due to the SO{3) rotational symmetry of the 3-space, 
we can decompose the metric perturbation into scalar, 
vector, and tensor modes. These modes are completely 
decoupled at linear order. We dehne the metric pertur¬ 
bation variables as 


<?oo = ~ (1 + 2a) , 

goi = a{t) {Si + diP) , 

g^j = a^{t) [5ij + 2'pSij + didjE 


The equation of motion for the tensor modes reads 

pk + 5Hjk+(^^ + ml^"fk = 0. ( 12 ) 

During inflation, the universe undergoes a superlumi¬ 
nal expansion with a nearly constant Hubble parame¬ 
ter. The vacuum fluctuations are stretched and frozen 
on super-horizon scales. At this stage, the graviton mass 
is 


rrpg ~ \p1{p^:/piY, because (/f, (13) 

where the subscript “f” is for “inflation”. If the gravi¬ 
ton mass were greater than the Hubble parameter iJ, 
the tensor modes would decay exponentially on super¬ 
horizon scales, and we would never see any signals today. 
To avoid the exponential suppression, the graviton mass 
must be much smaller than H during inflation, 


ml^XPliP./P.f^^Hl (14) 

As we shall see below, this condition can be easily satis¬ 
fied in our model. In passing, we note that as the scalar 
and vector modes are the same as in GR, our theory is 
free from the Higuchi ghost even for nig < 2H in the 
de-sitter space-time. 

Given the small but non-vanishing mass, the inflation¬ 
ary tensor spectrum is calculated by 


, _ 2h^ f k 

“ 7r2M| 


(15) 


with the tilt 


rit ~ —2e -|- 


2ml 


( 8 ) 


(16) 
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where e = —H/H^ is the slow-roll parameter. 

At the end of inflation, the slow-roll condition breaks 
down and the universe undergoes reheating. At the re¬ 
heating stage, the inflaton oscillates around the potential 
minimum, and gradually decays to radiation. The poten¬ 
tial is expanded around the minimum as 

+ ■■■ . (17) 

where M is the mass of the inflaton during reheating 
and the dots stand for higher order corrections which are 
irrelevant at low energy scale. Asymptotically for large 
Mt ^ 1, we have 


\/^Mt 


sin {Mt) exp ( --TMt ) , 


(18) 


where subscript “r” stands for “reheating”. The Ein¬ 
stein equations tell us 77 ~ ^ at this stage. To include 
the effect of decaying inflaton, we have simply added a 
decaying factor into the above solution, without 

specifying the detailed model of reheating. 

During the reheating stage, we have 


and thus (19) 

The equation of motion of gravitational waves (1121) be¬ 
comes a Mathieu-type equation, 


(Plk 

dx"^ 


-h 


2 d'jk , i-e 


-Tx 


dx 


■ sin^(x)7fc = 0, 


( 20 ) 


\ , 2 

where x = Mt and ^ ■ Noted that we have ne¬ 

glected the spatial gradient term since we are interested 
in the long wavelength modes. 

It is well known that the Mathieu equation has a very 
efficient and broad parametric resonance if ^ 1, 

i.e. the graviton mass must be much greater than Hub¬ 
ble parameter during reheating. Let us check whether 
this condition can be satisfied. Note that the Friedmann 
equation tells us MpH^ ^ M'^cj)^, where Hr is the Hubble 
parameter during reheating. Generally we expect that 
M^ ~ Hf due to the breaking of the slow-roll condition 
at the end of inflation. We thus get 


M^ 

Hi 


( 21 ) 


Demanding that the graviton mass be much greater than 
the Hubble parameter during reheating yields the condi¬ 
tion, 

M2 

Xj^ » 1. (22) 

Combinning conditions (fT4)) and (l2^ . we get 


M2 M2 



(23) 



FIG. 1: The resonant amplification of tensor modes during 
reheating. The horizontal axis is a: = Mt and the vertical 
axis is the relative amplitude of the tensor modes yfe/yj,*, 
where 7 *,* is the amplitude at the end of inflation. The pa¬ 
rameters are ^ = 10® and T = 0.05, with the initial condition 
d')k/dx\:r=l = 0 . 
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FIG. 2: The parameter dependence of the resonant amplifica¬ 
tion. The horizontal axis is F for the upper panels and ^ for 
the lower panels. The vertical axes is Logio{'yk/'yk*) evaluated 
at a: = 1000 . 


(j)- 

Thus ^4 ‘ 1, i.e. Mp is required for the self- 

consistency of the above inequality. Note that cj)^ -C Mp 
is also the condition of the validity of our effective field 
theory, which is automatically satisfied for many small 
field inflationary models. 

We have numerically solved Eq. (EOl)- The results are 
plotted in Figs. [Hand[2] The resonant amplification fac¬ 
tor depends on the value of ^ and the decay rate T. We 
can also read off the threshold for a significant resonant 
amplification is roughly ^ > 10^. The tensor modes stop 
growing in the large Mt limit due to the decay of the 
inflaton. For those long wavelength modes whose gradi¬ 
ent term is always negligible during reheating, the final 
power spectrum is still almost scale-invariant, as long as 
the graviton mass during inflation is small enough and 
thus the tensor tilt in Eq. (fTH]) is small. 

Noted that the kinetic term of inflaton field is canonical 
at leading order, which implies that in the massless limit 
TOg 0, scalar and tensor just simply decouple. On the 
other hand, it has been previously proven that this theory 
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smoothly reduces to GR in the massless limit due to the 
absence of vDVZ discontinuity 12|. After reheating, (f) —>• 
0, the graviton becomes massless and we recover GR. 

Vector perturbation To calculate the vector per¬ 
turbation, we adopt the unitary gauge, in which the fluc¬ 
tuations of 50(3) Stuckelberg scalar fields are fixed to 
be zero, i.e. <5(^9* = 0. The quadratic action of the vector 
perturbation reads (in momentum space) 


c(2) _ 

" 16 




FiFi - m^gF^Fi - 


4S,F, , 4SA 


(24) 


After integrating out Si, we get 


Sy^ - —Mp-ml 


a^k'^F.F,. 


(25) 


This clearly shows that the kinetic term for vector pertur¬ 
bation was canceled out. It is by no mean of an accident, 
because the kinetic term of vector modes are prohibited 
by internal symmetry —>■ -|- 5* . 

Scalar perturbation In the scalar sector, a, /3 and 
E are non-dynamical. After integrating them out, the 
quadratic action for the scalar perturbation in the uni¬ 
form (f) gauge (i.e. 64> = 0) reads 

Si« = A 4 " I ah (4 - , ( 26 ) 


This is exactly the same as the one in GR with a single 
scalar field. In this gauge, ijj is identical to the curvature 
perturbation on the comoving slicing, TZr., and the power 
spectrum is given by the same formula 21, ^ , 


Pn 


87r2eAf2' 


(27) 


Thus the tensor-to-scalar ratio we observe today is 

r = = i6e X A, (28) 

Pn 

where A is the resonant amplification factor of the tensor 
modes during reheating. For instance, with the parame¬ 
ters choice in Fig. [T] the factor A could be of the order 
of 10^^. The variation of the inflaton per e-fold is 



Thus during 60 e-folds, (j) traverses a distance Acj) ~ 
15Mpy^2r/A. Hence a sizeable tensor-to-scalar ratio is 
possible even for a sub-Planckian excursion. We conclude 
that the Lyth bound can be explictly evaded. 

Conclusion and discussion In this letter, we have 
minimally extended GR to a theory with a non-vanishing 
graviton mass term and proposed a mechanism to en¬ 
hance the primordial tensor perturbation from inflation 


in the the early universe. In our model, only the ten¬ 
sor perturbation is affected, while the scalar and vector 
perturbations remain the same as the ones in GR. The 
graviton mass is assumed to be proportional to the infla¬ 
ton during reheating, and hence its coherent oscillations 
give rise to a significant resonant amplification for all long 
wavelength modes on super-horizon scales. Then we have 
numerically studied the dependence of the amplification 
factor on the graviton mass and the decay rate of inflaton 
during reheating. We find that the Lyth bound can be 
explicitly evaded in our model. 


Our model contains three non-dynamical spacelike 
Stuckelberg fields which may formally become dy¬ 
namical if we include higher order derivative terms. Dur¬ 
ing inflation, however, the would-be new degrees of free¬ 
dom are supermassive and exponentially decay away. 
Thus we can safely integrate out these modes at low en¬ 
ergy scale. On the other hand, to screen these would- 
be new degrees at late time, if necessary, we can simply 
add a tiny but non-zero constant to the mass term in 
Eq. (ini). The current upper bound of the graviton mass 
nig is about I0“^°eV, from observation of Hulse-Taylor 
binary pulsar, PSR B I9I3-I-I6 23|. 


At the nonlinear perturbation level, we expect that the 
graviton mass term will introduce several new interaction 
terms. It will be interesting to study its possible imprints 
in, e.g. the non-Gaussianity of CMB anisotropies. As for 
models beyond the minimal model, a massive graviton 
generically induces non-trivial scalar and vector pertur¬ 
bations. We plan to study such possibilities and their 
possible observational effects in future. 


acknowledgments We would like to thank A. De Fe¬ 
lice, X. Gao, S. Mukohyama, T. Tanaka, A. Taruya, and 
N. Tsamis for useful discussions. 


[1] A. H. Guth, Phys. Rev. D 23 , 347 (1981). 

[2] Planck Collaboration, arXiv:1502.02114 [astro-ph.CO]. 

[3] P. Creminelli, D. L. Nacir, M. Simonovi, G. Trevisan and 
M. Zaldarriaga, arXiv:1502.01983 [astro-ph.CO]. 

[4] D. H. Lyth, Phys. Rev. Lett. 78, 1861 (1997). 

[5] M. Fierz, W. Pauli, Proc. Roy. Soc. Lond. A173, 211-232 
(1939). 

[6] C. de Rham, G. Gabadadze and A. J. Tolley, Phys. Rev. 
Lett. 106 , 231101 (2011) ;arXiv:1011.1232 [hep-th]]. 

[7] H. van Dam, M. J. G. Veltman, Nucl. Phys. B22, 397-411 
(1970). 

[8] V. I. Zakharov, JETP Lett. 12, 312 (1970). 

[9] D. G. Boulware, S. Deser, Phys. Rev. D6, 3368-3382 
(1972). 

[10] A. E. Gumrukcuoglu, C. Lin and S. Mukohyama, JCAP 
1203 , 006 (2012) arXiv:1111.4107 [hep-th]]. 

[11] A. De Felice, A. E. Gumrukcuoglu and S. Mukohyama, 
Phys. Rev. Lett. 109 , 171101 (2012). 

[12] S. L. Dubovsky, JHEP 0410, 076 (2004). 

[13] V. A. Rubakov, hep-th/0407104 















5 


[14] D. Comelli, F. Nesti and L. Pilo, JHEP 1307, 161 (2013) 
arXiv: 1305.0236 [hep-th]]. 

[15] b. Comelli, F. Nesti and L. Pilo, JCAP 1411, no. 11, 
018 (2014) arXiv:1407.4991 [hep-th]]. 

[16] C. Lin, Phys. Lett. B 727, 31 (2013) arXiv: 1305.2069 . 

[17] C. Lin, Phys. Lett. B 738, 386 (2014) \arXiv: 1307.2574 . 

[18] C. Lin and L. Z. Labun, arXiv:1501.07160 [hep-th]. 

[19] S. L. Dubovsky, P. G. Tinyakov and 1. 1. Tkachev, Phys. 


Rev. Lett. 94, 181102 (2005) hep-th/0411158 . 

[20] A. Higuchi, Nucl. Phys. B 282, 397 (1987). 

[21] M. Sasaki, Prog. Theor. Phys. 76, 1036 (1986). 

[22] V. F. Mukhanov, Sov. Phys. JETP 67, 1297 (1988) [Zh. 
Eksp. Teor. Fiz. 94N7, 1 (1988)]. 

[23] P. J. Sutton and L. S. Finn, Class. Quant. Grav. 19, 1355 
(2002) gr-qc/0112018 . 



